Caracterização de rearranjos cromossômicos citogeneticamente equilibrados associados a quadros clínicos

Este estudo teve como objetivos (a) identificar mecanismos pelos quais rearranjos cromossômicos citogeneticamente equilibrados possam estar associados de maneira causal a determinados quadros clínicos e (b) contribuir para a compreensão dos mecanismos de formação desses rearranjos. Para isso, foram estudados 45 rearranjos cromossômicos citogeneticamente equilibrados (29 translocações, 10 inversões e seis rearranjos complexos), detectados em pacientes que apresentavam malformações congênitas, comprometimento do desenvolvimento neuropsicomotor ou déficit intelectual. Foram 31 rearranjos cromossômicos esporádicos, três familiais que segregavam com o quadro clínico e mais 11 rearranjos cromossômicos herdados de genitores fenotipicamente normais. Inicialmente os pontos de quebra desses rearranjos foram mapeados por hibridação in situ fluorescente (FISH). A busca por microdeleções e duplicações genômicas foi realizada por a-CGH. A investigação dos pontos de quebra prosseguiu com a aplicação da técnica de Mate-Pair Sequencing (MPS), que permite localizar as quebras em segmentos de 100 pb - 1 kb, na maioria dos casos. Para obter os segmentos de junção das quebras no nível de pares de bases, os segmentos delimitados por MPS foram sequenciados pelo método de Sanger. A análise por aCGH revelou microdeleções ou microduplicações localizadas nos cromossomos rearranjados, em 12 dos 45 pacientes investigados (27%). A análise de 27 rearranjos por MPS permitiu a caracterização dos pontos de junção das quebras. MPS expandiu o número de pontos de quebra, detectados por análise do cariótipo ou aCGH, de 114 para 156 (em resolução < 2kb, na maioria dos casos). O número de pontos de quebra/rearranjo variou de 2 a 20. Os 156 pontos de quebra resultaram em 86 variantes estruturais equilibradas e outras 32 variantes não equilibradas. Perdas e ganhos de segmentos submiscroscópicos nos cromossomos rearranjados constituíram a principal causa ou, provavelmente, contribuíram para o quadro clínico de 12 dos 45 pacientes. Em cinco desses 12 rearranjos foram detectadas por MPS a interrupção de genes já relacionados à doença, ou provável alteração de sua região reguladora, contribundo para o quadro clínico. Em quatro dos 33 rearranjos não associados a perdas ou ganhos de segmentos, a análise por MPS revelou a interrupção de genes que já foram anteriormente relacionados a doenças, explicando-se, assim, as características clínicas dos portadores; outro rearranjo pode ter levando alteração da expressão gênica de gene sensível a dosagem e ao quadro clínico. Um rearranjo cromossômico familial, identificado na análise após bandamento G como uma translocação equilibrada, t(2;22)(p14;q12), segregava com quadro de atraso do desenvolvimento neuropsicomotor e dificuldade de aprendizado associados a dismorfismos. A combinação das análises por FISH, aCGH e MPS revelou que se tratava, na verdade, de rearranjo complexo entre os cromossomos 2, 5 e 22, incluindo 10 quebras. A segregação de diferentes desequilíbrios submicroscópicos em indivíduos afetados e clinicamente normais permitiu a compreensão da variabilidade clínica observada na família. Rearranjos equilibrados detectados em indivíduos afetados, mas herdados de genitores clinicamente normais, são, em geral, considerados como não tendo relação com o quadro clínico, apesar da possibilidade de desequilíbrios cromossômicos gerados por permuta desigual na meiose do genitor portador do rearranjo. Neste trabalho, a investigação de 11 desses rearranjos por aCGH não revelou perdas ou ganhos de segmentos nos cromossomos rearranjados. No entanto, a análise por aCGH da portadora de um desses rearranjos - inv(12)mat - revelou deleção de 8,7 Mb no cromossomo 8, como causa de seu fenótipo clínico. Essa deleção estava relacionada com outro rearranjo equilibrado também presente em sua mãe, independente da inversão. Para compreender os mecanismos de formação de rearranjos citogeneticamente equilibrados, investigamos os segmentos de junção no nível de pares de base. A análise por MPS que levou, na maioria dos casos, ao mapeamento dos pontos de quebras em segmentos <1kb permitiu o sequenciamento pelo método de Sanger de 51 segmentos de junções de 17 rearranjos. A ocorrência de blunt fusions ou inserções e deleções <10 pb, e a ausência de homologia ou a presença de micro homologia de 2 pb a 4 pb de extensão indicaram o mecanismo de junção de extremidades não homólogas (non-homologous end joinging; NHEJ), na maioria das 51 junções caracterizadas. As características de três dos quatro rearranjos mais complexos, com 17-20 quebras, indicaram sua formação pelo mecanismo de chromothripsis. Este estudo mostra a importância da análise genômica de variações de número de cópias por microarray, juntamente com o mapeamento dos pontos de quebra por MPS, para determinar a estrutura de rearranjos cromossômicos citogeneticamente equilibrados e seu impacto clínico. O mapeamento dos segmentos de junção por MPS, permitindo o sequenciamento pelo método de Sanger, foi essencial para a compreensão de mecanismos de formação desses rearranjos

Molekulare Analyse des Kolumnar-Gens beim Apfel (Malus x domestica)

Das Kolumnarwachstum beim Apfel (Malus x domestica) geht auf eine in den frühen 1960er Jahren entdeckte Zufallsmutation zurück. Die daraus resultierende Sprossmutante ist von großem wirtschaftlichem Interesse, da diese sehr kompakte Wuchsform unter anderem zu einer enormen Ertragssteigerung durch eine hohe Pflanzdichte der Bäume führt. Das Ziel der Arbeit ist die Entschlüsselung der molekularen Ursache dieser Mutation, die bisher weitgehend ungeklärt ist. Die Analyse wurde durch die Erstellung einer Referenzsequenz der Co-Zielregion einer kolumnaren Apfelsorte sowie durch die Konstruktion eng gekoppelter molekularer Marker realisiert. Durch die Konstruktion von genomischen Apfel-BAC-Bibliotheken mit mehrfacher Genomabdeckung und die Erstellung geeigneter Sonden wurde die Co-Region kloniert und deren Sequenz bestimmt. In Kombination zu dieser klassischen positionellen Klonierungsstrategie wurden genomische Illumina „mate pair“-Bibliotheken erstellt, sequenziert und bioinformatisch analysiert, um die genomische Region vollständig zu annotieren. Somit wurde eine vollständige genomische Referenz der Co-Region einer kolumnaren Apfelsorte erstellt, die die Grundlage für weitere Analysen bildet. Auf Basis dieser Referenz konnte die Co-Mutation in Form der Integration des LTR-Retrotransposons Gypsy-44 im kolumnaren Chromosom an Position 18,79 Mbp auf Chromosom 10 lokalisiert werden. Darüber hinaus konnten Transposon-basierende molekulare Marker erstellt werden, die eine verlässliche Genotypisierung von Apfelbäumen in Bezug auf das Kolumnarwachstum ermöglichen und dies unabhängig von der verwendeten Apfelsorte. Der genaue Wirkmechanismus von Gypsy-44, der zur Ausprägung dieses extremen Phänotyps führt, ist bislang unklar. Zusammenfassend lässt sich sagen, dass die molekulare Ursache für das kolumnare Wachstum aufgeklärt werden konnte und zudem die ersten molekularen Marker erstellt wurden, die eine sortenunabhängige Differenzierung zwischen kolumnaren und nicht kolumnaren Apfelbäumen ermöglichen.

Unlocking pathology archives for molecular genetic studies: a reliable method to generate probes for chromogenic and fluorescent in situ hybridization

Chromogenic (CISH) and fluorescent ( FISH) in situ hybridization have emerged as reliable techniques to identify amplifications and chromosomal translocations. CISH provides a spatial distribution of gene copy number changes in tumour tissue and allows a direct correlation between copy number changes and the morphological features of neoplastic cells. However, the limited number of commercially available gene probes has hindered the use of this technique. We have devised a protocol to generate probes for CISH that can be applied to formalin-fixed, paraffin-embedded tissue sections (FFPETS). Bacterial artificial chromosomes ( BACs) containing fragments of human DNA which map to specific genomic regions of interest are amplified with phi 29 polymerase and random primer labelled with biotin. The genomic location of these can be readily confirmed by BAC end pair sequencing and FISH mapping on normal lymphocyte metaphase spreads. To demonstrate the reliability of the probes generated with this protocol, four strategies were employed: (i) probes mapping to cyclin D1 (CCND1) were generated and their performance was compared with that of a commercially available probe for the same gene in a series of 10 FFPETS of breast cancer samples of which five harboured CCND1 amplification; (ii) probes targeting cyclin-dependent kinase 4 were used to validate an amplification identified by microarray-based comparative genomic hybridization (aCGH) in a pleomorphic adenoma; (iii) probes targeting fibroblast growth factor receptor 1 and CCND1 were used to validate amplifications mapping to these regions, as defined by aCGH, in an invasive lobular breast carcinoma with FISH and CISH; and (iv) gene-specific probes for ETV6 and NTRK3 were used to demonstrate the presence of t(12; 15)(p12; q25) translocation in a case of breast secretory carcinoma with dual colour FISH. In summary, this protocol enables the generation of probes mapping to any gene of interest that can be applied to FFPETS, allowing correlation of morphological features with gene copy number.

Deep RNA sequencing at single base-pair resolution reveals high complexity of the rice transcriptome

Understanding the dynamics of eukaryotic transcriptome is essential for studying the complexity of transcriptional regulation and its impact on phenotype. However, comprehensive studies of transcriptomes at single base resolution are rare, even for modern organisms, and lacking for rice. Here, we present the first transcriptome atlas for eight organs of cultivated rice. Using high-throughput paired-end RNA-seq, we unambiguously detected transcripts expressing at an extremely low level, as well as a substantial number of novel transcripts, exons, and untranslated regions. An analysis of alternative splicing in the rice transcriptome revealed that alternative cis-splicing occurred in similar to 33% of all rice genes. This is far more than previously reported. In addition, we also identified 234 putative chimeric transcripts that seem to be produced by trans-splicing, indicating that transcript fusion events are more common than expected. In-depth analysis revealed a multitude of fusion transcripts that might be by-products of alternative splicing. Validation and chimeric transcript structural analysis provided evidence that some of these transcripts are likely to be functional in the cell. Taken together, our data provide extensive evidence that transcriptional regulation in rice is vastly more complex than previously believed.

Genomic architecture of ecologically divergent body shape in a pair of sympatric crater lake cichlid fishes

Determining the genetic bases of adaptations and their roles in speciation is a prominent issue in evolutionary biology. Cichlid fish species flocks are a prime example of recent rapid radiations, often associated with adaptive phenotypic divergence from a common ancestor within a short period of time. In several radiations of freshwater fishes, divergence in ecomorphological traits - including body shape, colour, lips and jaws - is thought to underlie their ecological differentiation, specialization and, ultimately, speciation. The Midas cichlid species complex (Amphilophus spp.) of Nicaragua provides one of the few known examples of sympatric speciation where species have rapidly evolved different but parallel morphologies in young crater lakes. This study identified significant QTL for body shape using SNPs generated via ddRAD sequencing and geometric morphometric analyses of a cross between two ecologically and morphologically divergent, sympatric cichlid species endemic to crater Lake Apoyo: an elongated limnetic species (Amphilophus zaliosus) and a high-bodied benthic species (Amphilophus astorquii). A total of 453 genome-wide informative SNPs were identified in 240 F-2 hybrids. These markers were used to construct a genetic map in which 25 linkage groups were resolved. Seventy-two segregating SNPs were linked to 11 QTL. By annotating the two most highly supported QTL-linked genomic regions, genes that might contribute to divergence in body shape along the benthic-limnetic axis in Midas cichlid sympatric adaptive radiations were identified. These results suggest that few genomic regions of large effect contribute to early stage divergence in Midas cichlids.

Multiple parents crossover operators: A new approach removes the overlapping solutions for sequencing problems

Maintaining population diversity throughout generations of Genetic Algorithms (GAs) is key to avoid premature convergence. Redundant solutions is one cause for the decreasing population diversity. To prevent the negative effect of redundant solutions, we propose a framework that is based on the multi-parents crossover (MPX) operator embedded in GAs. Because MPX generates diversified chromosomes with good solution quality, when a pair of redundant solutions is found, we would generate a new offspring by using the MPX to replace the redundant chromosome. Three schemes of MPX will be examined and will be compared against some algorithms in literature when we solve the permutation flowshop scheduling problems, which is a strong NP-Hard sequencing problem. The results indicate that our approach significantly improves the solution quality. This study is useful for researchers who are trying to avoid premature convergence of evolutionary algorithms by solving the sequencing problems.

Extraction, PCR amplification, and sequencing of mitochondrial DNA from scent mark and feces in the giant panda

To expand the feasibility of applying simple, efficient, non-invasive DNA preparation methods using samples that can be obtained from giant pandas living in the wild, we investigated the use of scent markings and fecal samples. Giant panda-specific oligonucleotide primers were used to amplify a portion of the mitochondrial DNA control region as well as a portion of the mitochondrial DNA cytochrome b gene and tRNA(Thr) gene region. A 196 base pair (bp) fragment in the control region and a 449 bp fragment in the cytochrome b gene and tRNA(Thr) gene were successfully amplified. Sequencing of polymerase chain reaction (PCR) products demonstrated that the two fragments are giant panda sequences. Furthermore, under simulated field conditions we found that DNA can be extracted from fecal samples aged as long as 3 months. Our results suggest that the scent mark and fecal samples are simple, efficient, and easily prepared DNA sources. (C) 1998 Wiley-Liss, Inc.

A whole-genome massively parallel sequencing analysis of BRCA1 mutant oestrogen receptor-negative and -positive breast cancers

BRCA1 encodes a tumour suppressor protein that plays pivotal roles in homologous recombination (HR) DNA repair, cell-cycle checkpoints, and transcriptional regulation. BRCA1 germline mutations confer a high risk of early-onset breast and ovarian cancer. In more than 80% of cases, tumours arising in BRCA1 germline mutation carriers are oestrogen receptor (ER)-negative; however, up to 15% are ER-positive. It has been suggested that BRCA1 ER-positive breast cancers constitute sporadic cancers arising in the context of a BRCA1 germline mutation rather than being causally related to BRCA1 loss-of-function. Whole-genome massively parallel sequencing of ER-positive and ER-negative BRCA1 breast cancers, and their respective germline DNAs, was used to characterize the genetic landscape of BRCA1 cancers at base-pair resolution. Only BRCA1 germline mutations, somatic loss of the wild-type allele, and TP53 somatic mutations were recurrently found in the index cases. BRCA1 breast cancers displayed a mutational signature consistent with that caused by lack of HR DNA repair in both ER-positive and ER-negative cases. Sequencing analysis of independent cohorts of hereditary BRCA1 and sporadic non-BRCA1 breast cancers for the presence of recurrent pathogenic mutations and/or homozygous deletions found in the index cases revealed that DAPK3, TMEM135, KIAA1797, PDE4D, and GATA4 are potential additional drivers of breast cancers. This study demonstrates that BRCA1 pathogenic germline mutations coupled with somatic loss of the wild-type allele are not sufficient for hereditary breast cancers to display an ER-negative phenotype, and has led to the identification of three potential novel breast cancer genes (ie DAPK3, TMEM135, and GATA4).

Global methylation state at base-pair resolution of the Caulobacter genome throughout the cell cycle.

The Caulobacter DNA methyltransferase CcrM is one of five master cell-cycle regulators. CcrM is transiently present near the end of DNA replication when it rapidly methylates the adenine in hemimethylated GANTC sequences. The timing of transcription of two master regulator genes and two cell division genes is controlled by the methylation state of GANTC sites in their promoters. To explore the global extent of this regulatory mechanism, we determined the methylation state of the entire chromosome at every base pair at five time points in the cell cycle using single-molecule, real-time sequencing. The methylation state of 4,515 GANTC sites, preferentially positioned in intergenic regions, changed progressively from full to hemimethylation as the replication forks advanced. However, 27 GANTC sites remained unmethylated throughout the cell cycle, suggesting that these protected sites could participate in epigenetic regulatory functions. An analysis of the time of activation of every cell-cycle regulatory transcription start site, coupled to both the position of a GANTC site in their promoter regions and the time in the cell cycle when the GANTC site transitions from full to hemimethylation, allowed the identification of 59 genes as candidates for epigenetic regulation. In addition, we identified two previously unidentified N(6)-methyladenine motifs and showed that they maintained a constant methylation state throughout the cell cycle. The cognate methyltransferase was identified for one of these motifs as well as for one of two 5-methylcytosine motifs.

Cohabitation with a B16F10 melanoma-bearer cage mate influences behavior and dendritic cell phenotype in mice

This study evaluated the effects of cohabitation with a B16F10 melanoma-bearer cage mate on behavior and immune functions in mice. Five different experiments were conducted. In each of them, the female mice were divided into two groups: control and experimental. One mouse of each control pair was kept undisturbed and called ""companion of health partner"" (CHP). One mouse of each experimental pair was inoculated with B16FI0 cells and the other, the subject of this study, was called ""companion sick partner"" (CSP). On Day 20 of cohabitation, behavior and immune parameters from CHP and CSP mice were analyzed. In comparison to the CHP, the CSP mice: (1) presented an increased general locomotion in the open field and a decreased exploration time and number of entries in the plus-maze open arms; (2) had an enhanced expression of the CD80 costimulatory molecule on Iab(+)CD11c(+) spleen cells, but no differences were found on lymph nodes cells; (3) presented an altered differentiation of bone marrow cells in the presence of GM-CSF, IL-4, and LPS in vitro, resulting in a lower percentage of Iab(+)CD80(+) cells; (4) had a deficit in the establishment of a Delayed Type of Hypersensitivity to ovalbumin, which was associated to an in vitro proliferation of an IL-10-producing lymphocyte subpopulation after ovalbumin stimulation. Corticosterone levels detected on Day 20 of cohabitation were similar in CHP and CSP mice. It is shown here that DCs phenotype in mice is affected by conditions associated with behavioral alterations indicative of an anxiety-like state induced by the cohabitation with a tumor-bearer conspecific. This phenomenon occurred probably through a nondependent corticosterone mechanism. (C) 2009 Elsevier Inc. All rights reserved.

Isolation breeds naivety : island living robs Australian varanid lizards of toad-toxin immunity via four-base-pair mutation

Since their introduction to the toad-free Australian continent cane toads (Bufo marinus) have caused a dramatic increase in naïve varanid mortality when these large lizards attempt to feed on this toxic amphibian. In contrast Asian–African varanids, which have coevolved with toads, are resistant to toad toxin. Toad toxins, such as Bufalin target the H1-H2 domain of the α1 subunit of the sodium-potassium-ATPase enzyme. Sequencing of this domain revealed identical nucleotide sequences in four Asian as well as in three African varanids, and identical sequences in all 11 Australian varanids. However, compared to the Asian–African varanids, the Australian varanids showed four-base-pair substitutions, resulting in the alteration in three of the 12 amino acids representing the H1-H2 domain. The phenotypic effect of the substitutions was investigated in human embryonic kidney (HEK) 293 cells stably transfected with the Australian and the Asian–African H1-H2 domains. The transfections resulted in an approximate 3000-fold reduction in resistance to Bufalin in the Australian HEK293 cells compared to the Asian–African HEK293 cells, demonstrating the critical role of this minor mutation in providing Bufalin resistance. Our study hence presents a clear link between genotype and phenotype, a critical step in understanding the evolution of phenotypic diversity.

Next-generation sequencing reveals large connected networks of intra-host HCV variants

Background: Next-generation sequencing (NGS) allows for sampling numerous viral variants from infected patients. This provides a novel opportunity to represent and study the mutational landscape of Hepatitis C Virus (HCV) within a single host.Results: Intra-host variants of the HCV E1/E2 region were extensively sampled from 58 chronically infected patients. After NGS error correction, the average number of reads and variants obtained from each sample were 3202 and 464, respectively. The distance between each pair of variants was calculated and networks were created for each patient, where each node is a variant and two nodes are connected by a link if the nucleotide distance between them is 1. The work focused on large components having > 5% of all reads, which in average account for 93.7% of all reads found in a patient. The distance between any two variants calculated over the component correlated strongly with nucleotide distances (r = 0.9499; p = 0.0001), a better correlation than the one obtained with Neighbour-Joining trees (r = 0.7624; p = 0.0001). In each patient, components were well separated, with the average distance between (6.53%) being 10 times greater than within each component (0.68%). The ratio of nonsynonymous to synonymous changes was calculated and some patients (6.9%) showed a mixture of networks under strong negative and positive selection. All components were robust to in silico stochastic sampling; even after randomly removing 85% of all reads, the largest connected component in the new subsample still involved 82.4% of remaining nodes. In vitro sampling showed that 93.02% of components present in the original sample were also found in experimental replicas, with 81.6% of reads found in both. When syringe-sharing transmission events were simulated, 91.2% of all simulated transmission events seeded all components present in the source.Conclusions: Most intra-host variants are organized into distinct single-mutation components that are: well separated from each other, represent genetic distances between viral variants, robust to sampling, reproducible and likely seeded during transmission events. Facilitated by NGS, large components offer a novel evolutionary framework for genetic analysis of intra-host viral populations and understanding transmission, immune escape and drug resistance.

The genome sequence of the plant pathogen Xylella fastidiosa: The Xylella fastidiosa consortium of the organization for nucleotide sequencing and analysis, Sao Paulo, Brazil

Xylella fastidiosa is a fastidious, xylem-limited bacterium that causes a range of economically important plant diseases. Here we report the complete genome sequence of X. fastidiosa clone 9a5c, which causes citrus variegated chlorosis - a serious disease of orange trees. The genome comprises a 52.7% GC-rich 2,679,305-base-pair (bp) circular chromosome and 'two plasmids of 51,158 bp and 1,285 bp. We can assign putative functions to47% of the 2,904 predicted coding regions. Efficient metabolic functions are predicted, with sugars as the principal energy and carbon source, supporting existence in the nutrient-poor xylem sap. The mechanisms associated with pathogenicity and virulence involve toxins, antibiotics and ion sequestration systems, as well as bacterium-bacterium and bacterium-host interactions mediated by a range of proteins. Orthologues of some of these proteins have only been identified in animal and human pathogens; their presence in X. fastidiosa indicates that the molecular basis for bacterial pathogenicity is both conserved and independent of host. At least 83 genes are bacteriophage-derived and include virulence-associated genes from other bacteria, providing direct evidence of phage-mediated horizontal gene transfer.

Real-time quantitative broad-range PCR assay for detection of the 16S rRNA gene followed by sequencing for species identification

Here we determined the analytical sensitivities of broad-range real-time PCR-based assays employing one of three different genomic DNA extraction protocols in combination with one of three different primer pairs targeting the 16S rRNA gene to detect a panel of 22 bacterial species. DNA extraction protocol III, using lysozyme, lysostaphin, and proteinase K, followed by PCR with the primer pair Bak11W/Bak2, giving amplicons of 796 bp in length, showed the best overall sensitivity, detecting DNA of 82% of the strains investigated at concentrations of < or =10(2) CFU in water per reaction. DNA extraction protocols I and II, using less enzyme treatment, combined with other primer pairs giving shorter amplicons of 466 bp and 342 or 346 bp, respectively, were slightly more sensitive for the detection of gram-negative but less sensitive for the detection of gram-positive bacteria. The obstacle of detecting background DNA in blood samples spiked with bacteria was circumvented by introducing a broad-range hybridization probe, and this preserved the minimal detection limits observed in samples devoid of blood. Finally, sequencing of the amplicons generated using the primer pair Bak11W/Bak2 allowed species identification of the detected bacterial DNA. Thus, broad-spectrum PCR targeting the 16S rRNA gene in the quantitative real-time format can achieve an analytical sensitivity of 1 to 10 CFU per reaction in water, avoid detection of background DNA with the introduction of a broad-range probe, and generate amplicons that allow species identification of the detected bacterial DNA by sequencing. These prerequisites are important for its application to blood-containing patient samples.

Inheritance of female mating preference in a sympatric sibling species pair of Lake Victoria cichlids: implications for speciation

Female mate choice has often been proposed to play an important role in cases of rapid speciation, in particular in the explosively evolved haplochromine cichlid species flocks of the Great Lakes of East Africa. Little, if anything, is known in cichlid radiations about the heritability of female mating preferences. Entirely sympatric distribution, large ecological overlap and conspicuous differences in male nuptial coloration, and female preferences for these, make the sister species Pundamilia pundamilia and P. nyererei from Lake Victoria an ideally suited species pair to test assumptions on the genetics of mating preferences made in models of sympatric speciation. Female mate choice is necessary and sufficient to maintain reproductive isolation between these species, and it is perhaps not unlikely therefore, that female mate choice has been important during speciation. A prerequisite for this, which had remained untested in African cichlid fish, is that variation in female mating preferences is heritable. We investigated mating preferences of females of these sister species and their hybrids to test this assumption of most sympatric speciation models, and to further test the assumption of some models of sympatric speciation by sexual selection that female preference is a single-gene trait. We find that the differences in female mating preferences between the sister species are heritable, possibly with quite high heritabilities, and that few but probably more than one genetic loci contribute to this behavioural speciation trait with no apparent dominance. We discuss these results in the light of speciation models and the debate about the explosive radiation of cichlid fishes in Lake Victoria.